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We investigate pion source sizes through the Yano-Koonin-Podgoretskh 
(YKP) parametrization for the Hanbury-Brown Twiss (HBT) effect in 
the CERN-SPS 158 GeV/A central collisions. We calculate two-particle 
correlation functions numerically based on a (3+l)-dimensional rela- 
tivistic hydrodynamics with a first order phase transition and analyze 
the pair momentum dependence of the HBT radii extracted from the 
YKP parametrization in detail. We find that even in the case of a 
first order phase transition, expansion and the surface dominant freeze- 
out make the source in the hydro dynamical model opaque significantly. 
Consequently, the interpretation of the temporal radius parameter as 
the time duration becomes unavailable for the hydro dynamical model. 

Pion interferometry is one of the promising tools to investigate the hot and 
dense matter created in the relativistic heavy ion collisions It is well known 
as Hanbury-Brown Twiss (HBT) effect that we can obtain a size of the source 
from a two-particle correlation function because of the symmetry of the wave 
function of the emitted bosons. In the experimental analyses, the size of the 
source is obtained as a fitting parameter in the Gaussian fit to the correla- 
tion function and its physical meaning depends on the Gaussian fitting func- 
tion. For central collisions, Yano-Koonin-Podgoretskii (YKP) parametrization 
is used as one of the fitting functions. In the YKP parametrization, three spa- 
tial parameters are interpreted directly as the transverse (perpendicular to the 
collision axis), longitudinal (parallel to the collision axis) and temporal (emis- 
sion duration) extents, respectively 0. In the case of relativistic heavy ion 
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collisions, the reaction is highly dynamical. The collective expansion, which 
contains some informations of an equation of state of the matter, takes place 
naturally For example, if there exists a first order phase transition between 
quark-gluon plasma (QGP) and hadronic gas, enlarged time duration due to 
the existence of the mixed phase may be observed jj. However, the collective 
expansion makes the meaning of the size ambigious. Therefore, we investigate 
the sizes in the YKP parametrization based on a hydrodynamical model with 
a first order phase transition 

Assuming a completely chaotic source for simplicity, the two-particle cor- 
relation function is given as 

L- (ki, k 2 ) - 1 + {dN/d 3 ki){dN/d 3 k2 y W 

where ki and k 2 are momenta of a emitted pair and /(ki, k 2 ) is the interference 
term. 

As usual, we introduce the relative momentum = k^—k^ and the average 
momentum = [k^ + fc 2 )/2 that satisfy the on-shell condition q^K^ = 0. 
Assuming the cylindrical symmetry, we can take = (K°, Kt, 0, K£). Hence, 
only three components of relative momenta are independent. In the YKP 



parametrization, q±_ = Jq% + q^, q\\ = q z , and q° are taken as the independent 
components. Then the fitting function is given as 

Ctf,K") = l + \exv{-Rl(K»)ql-Rl(K»)[qf-(q ) 2 } 

- ~R 2 (K») + R 2 (K»)] [qX(K»)} 2 } , (2) 

where u M is a four- velocity which has only longitudinal component v(K^) in- 
terpreted as a expansion velocity and the fitting parameters are v(K^) and Ri 
(i = _L, ||,0). Because the three size parameters are invariant under longitu- 
dinal boosts, we may use a special frame called YKP frame, where v = 0, in 
the discussion of the meanings of the size parameters. Introducing the source 
function as 

s(x , m _ f U^)d^(x') U v {x>)k» 

b[X > K) -h (2vr)3 eMU P (x')kP/T)-l 6{X Xh {6) 
and the weighted average 

_ Jd*xA(x»)S(x»,K») 
{A[X }} ~ Jd*xS(x»,K») ' (4) 



the size parameters in the YKP frame are expressed as 



Ri(K») = (Ay)', (5) 
R 2 (Kn = (Atf - ^(xt) + ^[(Axf - (Ay)\ (6) 

P± P± 



JI {XZ) + J* 



R\(K») = (Azf - ^(x~z) + ^[(Ax) 2 - (Ay)% (7) 



where x — x — (x), Ax = yj (x 2 ) — (x) 2 , j3± = K T /K° and (3\\ = K L /K°. From 
the expressions above, we can interprete the size parameters as the extents of 
the source if the last two terms in Eqs. @ and ([F|) are sufficiently small. It 
has been already shown that the interpretation holds within a class of thermal 
source model 0, and the case of a hydrodynamical model is the point of the 
present paper. 

Figure 1 shows dependence of the pion HBT radii. Close circles stand 
for the HBT radii extracted from the pion correlation function (|l|) through 
the fitting. Solid lines are space-time extensions calculated from Eqs. (§j)-(0), 
which are expected to agree with the HBT radii when only thermal pions are 
considered. Dotted lines stand for the "source sizes" (Az for R» and At for 
Rq). As far as R± and R\\ are concerned, those quantities seem to be consistent 
and well reproduce the NA49 experimental results || (open circles). As for 
Rq, our results are a few fm smaller than the experimental results in spite 
of the existence of the first order phase transition. Furthermore, the time 
durations At do not agree with the Ro at large Kt, that means the failure of 
the interpretation of Rq as the time duration. As explained below, this failure 
comes from an opaque property of the source. 
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Fig. 1 Kt dependence of YKP radii. Open and closed circles show the ex- 
perimental data and HBT radii, respectively. Solid lines stand for space-time 
extensions Dotted lines stand for source sizes. (Az for Ru and At for 

Ro.) 



The "opaque" source || emits particles dominantly from the surface and 
the "transparent" source emits from whole region. Figure 2 shows the source 
function ([|) projected onto the x — y plane, St(x,u) = * J dzdtS(x^ , K^) be- 
ing normalized as / dxdyST{x,y) = 1. Note that the average momenta of the 
emitted pions is x direction in the figure. The left figure shows clearly surface 
emission; pions are emitted from the crescent region and the source is thin- 
ner in the x direction than in the y direction. This is typical property of the 
opaque source and should be appeared in the third and fourth term in Eq. (^) 
and the factor (Ax) 2 — {Ay) 2 works as a measure of source opacity. On the 
other hand, if we neglected the transverse flow by putting U r = by hand, the 
source function becomes as shown in the right figure of Fig. 2. In this case, 
the source function is proportional to the space-time volume of freeze-out hy- 
persurface. As a consequence of restoration of the azimuthal symmetry which 
is shown clearly in the figure, the above measure of source opacity vanishes. 
Although the particle emission takes place almostly from the thin surface, the 
source is not opaque in this sence. When the transverse flow exists, the source 
function is deformed by thermal Boltzmann factor exp(Kj< smh. Yt cos (f)/Tf) 
with Yt almost proportional to r. Consequently, the number of emitted parti- 
cles increases in the region where x > (i.e., cos0 > 0) and decreases in the 
region where x < (cos0 < 0). This flow effect deforms the surface dominant 
volume distribution (right figure) to the crescent shape (left figure). 
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Fig. 2 The source functions as functions of transverse coordinates x and y for 
K T = 450 MeV and Y wn = 4.15 (Y n7T = § In guffi )- The left fi g ure denotes 
the source function with the transverse flow. The right figure denotes the 
source function without transverse flow. 

In summary, we analyze dependence of source parameters of the YKP 



parametrization based on the relativistic hydrodynamics for the CERN-SPS 
158 GeV/A Pb+Pb collisions. We obtain the results almostly consistent with 
the experiment. However, the temporal source parameter, Rq, shows the be- 
havior different from the experiment. The source opacity makes the interpre- 
tation of Ro as the time duration hard. We found that the source opacity was 
caused by the transverse flow and the characteristics of the freeze-out hyper- 
surface, the surface dominant freeze-out. The deviations of our results from 
the experiment would be improved by including the resonance decay and other 
effects. 
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